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of associative ionization have been measured using an original method based on relative 
measurements of NaJ and Na+ ion signals. 
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1. Introduction 

Collisional ionization of atoms may often be associated with highly excited quasi- 
molecular systems, which are formed temporarily in course of the collision. Such quasi- 
molecules may, at certain internuclear distances, couple to the ionization continuum and 
emit an electron. Depending on the total excitation energy of the system and relative 
positioning of the covalent and ionic molecular potential curves, the quasi-molecule may 
end up as a stable molecular ion (associative ionization), or ionize and dissociate into 
an atomic ion and a neutral atom (Penning process). Conventional theory ^ 12 E] 
treats such processes deterministically in terms of interaction of a few covalent bound 
states with the ionization continuum. In the reality, however, the incoming covalent 
state interacts with a large number of other molecular Rydberg states. Moreover, since 
the density of Rydberg states is very high, it may be impossible to single out separate 
interactions between two Rydberg states at a time. Obviously, under such conditions it 
is not practicable to apply the exact quantum mechanical description. This drawback 
of the deterministic description of ionization dynamics can be overcome by involving a 
new alternative approach, the so-called stochastic ionization model 

The stochastic model describes the collisions with Rydberg atoms as a chaotic 
migration of highly excited electron accompanying the motion of both colliding nuclei. 
The electron "diffuses" through the dense energy spectrum of the highly excited quasi- 
molecule, eventually ending up in the ionization continuum. The stochastic theory 
has already yielded encouraging results by properly describing some of the earlier 
experimental results on associative ionization involving Rydberg atoms IHl E] . 

The aim of the present work was to test experimentally the ranges of validity of 
the stochastic theory developed in [H El 13 IE] • For this purpose, we studied the 
ionizing collisions of sodium Rydberg atoms Na(nZ/) with atoms in the 3S ground state 
Na(3S). In what follows we describe the measurements of ionization signals occurring 
after excitation of sodium atoms to nS and nD Rydberg states with n=8-20 in a single 
atomic beam at a source temperature of 635 K. To the best of our knowledge, there are 
no data in the literature on ionization of sodium atoms in nS and nD states for this 
range of principal quantum numbers. In particular, the processes of interest are: 

Na(?T,L) + Na(3S) Na^ -|- e Associative Ionization (AI); 

Na(nL) + Na(3S) Na+ + Na(3S) + e Penning Ionization (PI); 

where Na"*" is an atomic ion, and NaJ is a molecular ion. 

The earlier experiments on collisional ionization of Na Rydberg atoms are reported 
in inmUmillEllinillllllSl- The measurements of absolute rate constants have been done 
in P for nP states {n = 5 — 15) in a crossed beam experiment at source temperatures 
of 600 K and pulsed laser excitation. The rate constants for nP states {n = 5 — 21) 
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Figure 1. Experimental data on the associative ionization rate constants of nS, nP, 
and nL{L > 1) Rydberg states of sodium (from refs. pill2|). 



in a single atomic beam at 700 K source temperature and continuous excitation were 
reported in fOlE]- The measurements for nS, nP, and nL{L > 1) states {n = 17 — 27) 
in a single Na beam at 1000 K and pulsed laser excitation have been made in [T2]. The 
use of pulsed excitation in jHl [T2j allowed the authors to distinguish between AI and PI 
channels and to determine the AI rate constants alone, while only the total (AI+PI) 
rates could be measured in ^01 E]- The compilation of the AI rates from jHl is 
shown in figure 1. It exhibits a maximum in the AI rate constant for Na(nP) states 
near n=12. Similar dependences on the principal quantum number were observed in 
experiments with other alkali Rydberg atoms pL^J. Comparison of these data with the 
theory of [11121 E] has shown a remarkable disagreement of the absolute values and shapes 
of experimental and theoretical n dependences. Therefore new experiments and further 
development of the theory are necessary. 

One can see from figure 1 that no experimental data for sodium nS and nD states 
with n < 17 are available. Meanwhile, the range of n = (8 — 20) is of particular interest, 
since this is where the maximum in the ionization rate is expected. In addition, S and 
D states have different quantum defects and orbits of electron motion from those of P 
states, and the n dependences for S and D states may be considerably different. 

One of the main difficulties in measurements of AI rates for nS and nD states at 
low n is related with their short lifetimes. According to the calculations of Theodosiou 
[T^ . radiative lifetimes of the lOS and lOD states are about 0.9 fis, while the lifetime 
of the lOP state is about 8 /is, which is by an order of magnitude larger. Since one 
must usually provide a sufficiently long time (1-2 /is) for the collisional interaction of 
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Rydberg and ground-state atoms in experiments, care must be taken to renormalize the 
ionization signals and thus to account for the short hfetimes. 

Experimental data on effective lifetimes of most of the Rydberg states are lacking, 
therefore they must be accurately calculated. Such calculations should take into account 
the fact that lifetimes are affected by blackbody radiation (BBR), which is always present 
in the experiments. It induces transitions between neighboring Rydberg states and 
shortens the effective lifetimes ^7] . Moreover, BBR may photoionize Rydberg atoms and 
contribute Na'*' ions to the measured collisional signals. Our numerical calculations have 
shown (see summary Table 1 in Section 4 below) that the rate of direct photoionization 
by BBR is as high as 1147 s~^ for the 17D state at the ambient temperature T=300 
K, so that the measured Na"*" signals in room-temperature experiments are rather due 
to the BBR photoionization than due to the PI. This effect was studied experimentally 
earlier in I14 | ITH]. while numerical calculations of direct BBR photoionization rates 
were performed in ^Hj . Experimental measurements of BBR photoionization rates 
for various n were reported only for barium Rydberg atoms |2()j . 

In order to decrease the BBR in the experiments with Rydberg atoms, special 
thermal shields cooled by vacuum cryostats are often used. However, this method 
does not allow a precise determination of the effective temperature and density of 
blackbody photons. Even with liquid nitrogen or helium cooled shielding, the blackbody 
photons can penetrate into the reaction zone from the hot atomic beam source (oven) 
or through openings in the detection system. Although our setup was equipped with a 
nitrogen cryostat, we decided to work without cooling for the sake of exact knowledge 
of temperature and frequency distribution of BBR photons. Under such conditions, 
reliable measurements of PI rates were impossible because the expected Na"*" signals 
due to PI were much smaller than those from BBR photoionization. 

At the same time, the BBR photoionization signals may be calculated or measured 
and serve as a reference for measurements of the AI rates. In this case, however, the 
mixing of neighboring Rydberg states by BBR must also be taken into account, because 
it may significantly redistribute the populations of excited states during the interaction 
time pn 122] ■ Note that the AI rates can still be extracted directly from the measured 
NaJ signals, since the BBR photoionization leads only to the formation of Na"*" ions and 
does not affect the Na^ signals. 

2. Experiment 

Experimental setup was similar to that used for microwave spectroscopy of Na Rydberg 
atoms [23]. Experiments were performed with a single effusive Na atomic beam in a 
vacuum chamber at background pressure of 5-10"'' Torr (figure 2). The temperature 
of the Na oven was stabilized at 635 K. The atomic beam was formed by expansion of 
sodium vapor through a 2 mm dia opening in the oven at a distance of 9 cm from the 
reaction zone. CoUimation of the beam was achieved by a metallic shield with 1.5 mm 
dia aperture located 4 cm upstream from the excitation volume. The effective diameter 
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Figure 2. Experimental arrangement of the atomic and laser beams, and the ion 
detection system. 

of the atomic beam in this volume was about 4 mm. 

The metalhc shield screened the most of thermal radiation of the oven. This shield 
and other parts of the vacuum chamber surrounding the reaction zone were kept at 
room temperature, ensuring an isotropic 300 K environment. The estimated flux of 
"hot" BBR photons originating from the beam source and penetrating through the 
collimating aperture contributed to less than 1% of the total BBR ionization rate in the 
reaction zone, and it was therefore disregarded in the data analyses. 

Sodium nS and nD Rydberg states were excited using the two-step scheme 
3Si/2 ^3P3/2 — >nS, nD by radiation of two tunable lasers pulsed at 5 kHz repetition 
rate. In the first step, 50 ns pulses from a Rhodamine 6G dye-laser with linewidth of 50 
GHz were used. They saturated the 3Si/2 ^3P3/2 transition at 589 nm. The resonance 
fluorescence on this transition was detected by a photomultiplier to monitor the relative 
changes in density of the atomic beam. In the second step, the second harmonic of a 
Ti-sapphire laser was used. It yielded 50 ns pulses with 10 GHz linewidth, tunable in the 
400-455 nm range. When not focused, this radiation did not saturate the 3P3/2 -^nS,nD 
transitions. The two laser beams were crossed at a right angle in the reaction zone, both 
of them crossing the atomic beam at an angle of 45°. Laser beam profiles were spatially 
limited by 2 mm dia apertures at the entrance windows of the vacuum chamber. Such 
configuration ensured a sufficiently small excitation volume of 2 mm size in the central 
part of the atomic beam, where the spatial variation of atom density was insignificant 
(<20%). 

The ion detection system shown in figure 2 used a channeltron multiplier VEU — 6. 
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The atomic beam passed between two stainless-steel plates with diameter of 70 mm 
each, spaced by 10 mm. The plates formed a homogeneous electric field to guide the 
ions from the reaction zone through the entrance window of the channeltron. The 
extraction electric field pulses of 100 V/cm amplitude and 250 ns duration were applied 
to the upper plate. The lower plate was grounded and had a 6 mm dia opening covered 
by a mesh with the transmittance of 70%. The ions that passed through the mesh were 
accelerated by the electric field of the channeltron to energies of about 3.5 keV. This 
energy was sufficient to detect the ions with probability close to unity. 

A special microwave-spectroscopy experiment, similar to that described in |23j . 
has proven that the electric field from the channeltron very weakly penetrated into the 
interaction region through the fine mesh of the lower ionization plate, as well as other 
stray electric fields. The measured Stark shift of the microwave transition 37S^37P 
was about 4 MHz, corresponding to the <0.3 V/cm stray electric field. This field had 
no effect on our measurements or on Rydberg states under study. 

Single ion output pulses of the channeltron were amplified, discriminated, and 
registered by two gated counters. The measurements were performed in the pulse 
counting regime, keeping the frequencies of the detected ion signals much lower (0.2- 
1 kHz) than the 5 kHz repetition rate of laser pulses. On average, less than one ion per 
laser shot was detected. The measured frequencies of ion signals were determined as the 
total number of ions detected during the measurement time of 10 s, i.e., signals were 
counted for 50000 laser pulses. In order to ensure the single-ion counting regime, the 
intensity of the laser driving the second excitation step was attenuated with calibrated 
neutral density filters by a factor of 10 to 100. 

The present study concentrated on measurements of relative dependence of 
ionization rate constants on the principal quantum number n of Rydberg states, and 
it did not require precise knowledge of absolute values of the number density n^s of 
the ground-state Na(3S) atoms in the beam. It was therefore calculated from the 
geometry of the effusive beam and the Nesmeyanov's formula [25 relating temperature 
and pressure of saturated Na vapor in the beam source. This formula coincides to 
within better than 50% with other available data, for example, [23]. Monitoring of 
the fluorescence on the saturated resonance transition showed that the atomic number 
density was almost constant during the experiments. The temperature of the oven was 
measured by a calibrated thermocouple with an accuracy of ±2 K. Taking into account 
all the uncertainties, we estimate the absolute number density of atoms in the reaction 
zone as n35'=(2±0.7) xlO^° cm^^ at the oven temperature of T=(635±2) K. 

The time sequence of excitation and detection pulses is illustrated in flgure 3. 
Two identical electric fleld pulses with the 100 V/cm amplitude and 250 ns duration 
(flgure 3(b)) were applied to the repelling plate after each laser excitation pulse 
(flgure 3(a)). 

The flrst electric pulse was applied immediately after the laser pulse to remove the 
atomic and molecular ions produced during the laser pulse. These ions were always 
observed in our experiments, even at non-resonant excitation by the flrst-step (yellow) 
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Figure 3. Timing diagram of signals: (a) laser excitation pulse; (b) 100 V/cm electric 
field pulses for ion extraction; (c) ion signals; (d) counter gates for measurement of 
Na^ and Na^ signals. 



laser only. Their origin was investigated in many earlier studies, which identified them 
to be due to mutiphoton ionization of Na2 dimers by yellow laser [211 |22I (dimers are 
always present in an effusive beam at an amount of ~(0.1-1)%), due to photoassociative 
ionization of Na(3P) atoms by yellow laser [23 1211 EOl; and due to photoionization of 
Rydberg Na.{nL) atoms by the resonant radiation of yellow and blue lasers. These ions 
represent an undesirable background to the collisional-ionization signals from Rydberg 
atoms, and therefore the first electric pulse was required to clean up the reaction zone 
from any ions created during the laser excitation. 

The second electric pulse extracted to the channeltron those ions, which were 
produced in the time interval between ti=0.3 /xs and ^2=2.1 /xs after the laser excitation 
pulse. These ions corresponded to the collisional and BBR photoionization of Rydberg 
atoms. In the mass spectrum detected by the channeltron, the signals of the atomic 
Na^ and the molecular Na^ ions were separated by 0.6 fis and thus well resolved 
(figure 3(c)). Gated pulse counters registered the signals of the atomic and molecular 
ions independently (figure 3(d)). 
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3. Determination of collisional rate constants 
3.1. Rate equations 

The registered Na"*" and Na^ ion signals result from ionization occurring in the reaction 
zone during the time interval {t2-ti)=1.8 /xs between the two extraction electric field 
pulses. This time is comparable with the lifetimes of Rydberg states; therefore, time 
evolution of ionization processes must be analyzed. The main processes leading to the 
formation of Na+ ions are Penning ionization and photoionization by BBR. The Na^ 
ions can be created only in the associative ionization. Contribution of the collisions with 
background gases can be safely disregarded. We have verified experimentally that the 
variation of background pressure within the range of 5-10-^ < P < 1-10'^ Torr did not 
affect the measured Na"*" and NaJ signals by more than 5%. Under such conditions, the 
rate equations describing the evolution of the numbers of Na"*" and Na^ ions following 
the laser excitation at time t=0 are 

= kpjNnL (t) ^35 + WsBRNnL (t) ; 

(1) 



dt 

Here, Nniit) is the time-dependent number of Na(nL) Rydberg atoms in the reaction 
zone, Ri2-10^° cm~^ is the number density of ground state atoms, kAi and kpi are 
the rate constants of associative and Penning ionization in Na(niy)-|-Na(3S) collisions 
respectively. 

The function NnL{t) evolves in time as 

NnL (t) ^ NnL (0) CXp [-t/nff] . (2) 

The initial number of Rydberg atoms, A'^„l(0), created during the laser excitation, 
can be written as 

NnL (0) = NspW (3P3/2 nL) , (3) 

where N^p is the average number of atoms in the 3P3/2 state during the yellow-laser 
pulse, and W (3P3/2 — > nL) is the probability of excitation of the Na(3P3/2) atoms to 
the nL state by a single blue-laser shot. 

The effective lifetime Tg// describing the decay of Rydberg states in equation Q 
is determined by the spontaneous lifetime and the rate of other processes depleting the 
laser excited Rydberg state. These include BBR induced transitions between Rydberg 
states, BBR induced photoionization, and collisional quenching. The results of our 
numeric calculations of the rates of these processes will be given in Section 4 and 
summarized in Table 1. 

The depletion of Rydberg states with n=8-20 by collisional ionization is negligible 
at atom densities used in the present experiment. We find from figure 1 that the rate of 
associative ionization, kAjn^s, does not exceed 50 s~^ and is therefore much smaller than 
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the spontaneous decay rates, which range from 10^ to 10^ s~^ for the studied Rydberg 
states. The rate of PI, kpjn^s, is expected to be below 10 s~^ for n~20, and close to zero 
for lower n. Comparing the PI rate with the direct BBR photoionization rate W bbr 
(Table 1), one can see that Na"*" ions are produced mainly in the BBR photoionization. 
As will be shown below, this seemingly undesirable background ionization process can 
be favorably exploited for the determination of absolute AI rate constants. 

With the above considerations in mind, the solution of equations (Q) can be written 

as 

Na+ = Ar,i(0) W BBR teff 

Na+ = NnL (0) kAin^s teff ^ ' 

where teff is the effective time of interaction that takes into account the short radiative 
lifetimes of Rydberg states: 

= ^eff [exp {-t J Teff) - exp {-t J Teff)] . (5) 

Here, ti and ^2 are the time moments marking the end of the first and the beginning of 
the second extraction electric field pulses, respectively (see figure 3(b)). 



3.2. Method of measurements 

Equations(jlI) can be used for the direct measurement of kAi and W bbr values, provided 
that Nnhif^) is known. The only reliable method to measure Nniif^) is the Selective 
Field Ionization (SFI) technique jSI]. Rydberg atoms ionize in an electric field with 
probability close to unity if the field strength has reached the critical value Ee- The 
latter depends strongly on the effective quantum number neff={n-6L), where 5^ is the 
quantum defect (1.348 for nS states and 0.015 for nD states): 

Ee ^ 3.2 • lO^n^ff (V/cm) (6) 

Unfortunately, it is difficult to apply the SFI method for Rydberg states with low n, since 
they require too strong electric field (~30 kV/cm for n~10). In the present study we 
were interested mainly in relative measurements of kAi for various n. Therefore we could 
use a normalization procedure for NnL{0) based on numerically calculated excitation 
probabilities W (3P3/2 — * nL). Since the 3Si/2 ^3P3/2 transition was saturated, A'^„l(0) 
was dependent only on the respective transition moments and power of the blue laser. 
In the absence of saturation at the second excitation step (this was the case in our 
experiments), the probability of excitation of Rydberg states from the 3P3/2 state can 
be written as 

W (3P3/2 ^ nL) =CL-Pb-R^ (3P3/2 ^ nL) , (7) 

where Ph is the power of the blue laser, R (3P3/2 nL) is the radial part of the 
transition dipole moment, and C l is a, normalization constant, which depends on L 
and is proportional to the square of angular part of the matrix element. The value of 
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W (3P3/2 — ^ nVj falls as n~jj for high Rydberg states, but for the states with n~10 this 
scaling law does not work quite well. This was established in our numeric calculations 
of R (3P3/2 nL) for the 3P3/2 ^nS, nD transitions, therefore we used the calculated 
values instead of the above scaling law for the evaluation of the experimental data. 

In order to compare the absolute signals due to BBR and collisional ionization of 
S states with those of D states, it is necessary to know also the ratio CojCs- The 
analysis of angular parts of the transition matrix elements, taking into account the 
hyperfine structure, showed that for excitation with linearly polarized light in the first 
and the second excitation steps, the ratio Co/Cs may vary from approximately 1.6 
(if there is no collisional, radiative, or magnetic field mixing of the magnetic sublevels) 
to 2 (if the level mixing is complete). For excitation by non-polarized light, the ratio 
always equals to 2 regardless the degree of level mixing. Finally, we find that the ratio 
W (3P3/2 — > nD) /W (3P3/2 nS) may vary between the values of 3.5 and 5. 

In principle, one could normalize the ion signals measured for different nL states 
using the calculated probabilities W (3P3/2 nL) and measuring only the power Pf, of 
the blue laser radiation in equation ((Tj). However, the applicability of such normalization 
may be complicated by technical imperfections of the blue laser. Since the linewidth 
of this laser (10 GHz) was much larger than the widths of the absorption profiles at 
the second excitation step (~500 MHz Doppler broadening), variations of the spectral 
density of laser radiation could affect the probability of excitation even if P;, would 
be kept constant. Therefore we had to verify experimentally the applicability of 
normalization by equation (jZj). As discussed above, the only reliable way to measure the 
number of Rydberg atoms was to apply the SFI technique. For this purpose, we built 
a high-voltage generator yielding pulses with rise time of 1 /is and amplitude of up to 
8 kV. This allowed us to field-ionize Rydberg states with n>17. The SFI signals were 
detected at a 1 /xs delay with respect to the laser pulse, i.e., the measured SFI signal 
was: 

~ A^nL(0)exp(-l /is/re//). (8) 

Equation® was used to derive NniiS^) from the measured SFI signals and the 
calculated values Te// of Table 1. Figure 4 shows the measured NnhiS^) dependences on 
the principal quantum number n for nS and nD states. These data are normalized on 
Pfe, because it varied as blue laser frequency was tuned to resonance with different nL 
states. The solid curves are approximations made using equation ((7j). It is seen that 
experimental points have noticeable deviations from theory although the general trend 
is correct. These deviations may be explained by the variations of spectral density of 
blue laser light. We concluded that equation ((Tj) can be used for the normalization of 
^raL(O), but at price of limited accuracy. We also find from figure 4 that on average the 
ratio iy(3P3/2 — >-nD)/iy(3P3/2 — >nS) was close to 3.5. Hence, no considerable mixing 
of the magnetic sublevels took place during the laser excitation, and the ratio Co/ Cs 
was close to 1.6. 
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Figure 4. Relative probabilities of laser excitation of sodium Rydberg states: (a) nS 
states; (b) nD states. Open circles and squares - experiment, solid curves - theory. 



Since SFI technique was not applicable for the direct determination of NnL{^) for 
n=8-17, the use of equation (0) was somewhat inadequate in our experiment. Therefore 
ways had to be found to eliminate this value from the determination of rate constants. 
This can be done by measuring the ratio R of atomic and molecular signals derived from 
equations (jH): 



R 



Na^ kAin^s 



Na^ 



Wbbr 



(9) 



This ratio does not depend on the values of NnL^O), r^jf and t^jf . Thus, the rate 
constant of the AI process can be directly obtained from the measured ratio of the Na^ 
and Na^ signals: 



k 



Nat W, 



Ar = p,-^. (10) 
Na^ ^35 

The BBR ionization rates Wbbr have thus become one of the key values necessary 
for the determination of the AI rate constants. The accuracy with which the Wbbr 
values are known directly affect the accuracy of the experimental kAi values obtained 
using equation (fTIUl . Therefore, the process of BBR ionization must be examined in all 
details. 



4. BBR ionization rates 



4.1. Direct BBR photoionization and effective lifetimes 

First, we have numerically calculated the effective lifetimes Tg// and the rates Wbbr 
of direct BBR photoionization of nS and nD Rydberg states of Na with n=8-20 at the 
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room temperature of 300 K. The BBR photoionization rates were calculated using the 
formula (here and below the atomic units are used) 



Wbbr = c j p^du, (11) 

where c is the speed of light, LOnL = 1/ (2^^e//) threshold frequency of the nL 

Rydberg state, and is the photoionization cross-section at frequency cu. The volume 
density p^^ of BBR photons at temperature T is given by the Planck distribution: 

= 7r2c3 [e-/(*^^) - 1] ' ^^^^ 
where k is the Boltzmann constant. For isotropic and non-polarized thermal radiation 
field, is determined by the radial matrix elements R {nL —^E,L±1) of dipole 
transitions from discrete nL Rydberg states to continuum states with L' = (L ± 1) 
and photoelectron energy E: 



^"^slilTTT 5Z L^a.li'{nL^EL±l), (13) 

^ ' L'=L±l 

where L^ax is the largest of L and L' . The values of R [uL —>■ E,L ±1) were calculated 
numerically using the semi-classical formulae of Dyachkov and Pankratov 32 . Their 
method gives orthogonal and normalized continuum wavefunctions, and the calculated 
photoionization cross-sections were found to be in good agreement with the more 
sophisticated quantum-mechanical calculations for nS states by Aymar [33] • Other 
theoretical methods developed in [SH ES] give close results. High-precision values of 
quantum defects of Dyubko et al jHE] were used in our calculations. 

In order to determine the effective lifetimes of uL states, the radial matrix elements 
R {nL — i> n'L') of all dipole transitions to other n'L' states with L' = L ± 1) were also 
calculated using the quasi-classical formulae of [32] ■ The rate of a specific BBR-induced 
transition between the states nL and n'L' are given by: 

W {nL ^ n'L') = A^I^R^ ^nL n'L') (14) 

^ ^ 3c3 2L + 1 ^ ^ e'^-'/('=^) - 1 ^ ^ 

where Unn' = l/(2n^) — l/(2n'^) is the frequency of the transition. The total decay 
rates were obtained by summing the probabilities of all spontaneous and BBR-induced 
transitions. The results of numerical calculations of Wbbr and effective lifetimes at 
T=300K are presented in Table 1. 



4-2. BBR-induced mixing of Rydberg states 

Earlier studies of BBR-induced processes have revealed that BBR causes not only 
direct photoionization to the continuum ^Hl, but also induces transitions in the discrete 
spectrum |211 122] • In our case, several discrete nP and nF levels are populated from 
the initially excited nS and nD states during the interaction time. For example, after 
laser excitation of the 16S state, the BBR-induced transitions populate the neighboring 
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Table 1. Effective lifetimes Tefj{fis) and BBR induced transition rates of nS and 
riD states of sodium for the surrounding environment at T=300 K. Wbbr - direct 
photoionization rates of laser excited Rydberg states; Wb. - ionization rates by the 
second extraction field pulse of states with n' > 40, which are populated in BBR 
induced transitions from the initial laser excited levels; Wmix - rate of BBR induced 
ionization of n'L' levels populated by spontaneous decay and BBR-induced nL mixing; 
^BBR ' total theoretical BBR ionization rate; W^^^ - experimental BBR induced 
ionization rate. 



State 




Wbbr{s-^) 




Wrnix{s 1) 




tttCXV / 1 \ 


8S 


0.395 


0,059 


0.022 


0.404 


0.485 




9S 


0.581 


0.823 


0.316 


3.62 


4.76 


6±2 


lOS 


0.812 


4.555 


1.78 


13.7 


20.0 


25±2 


lis 


1.095 


14.48 


5.76 


30.9 


51.2 


47±3 


12S 


1.431 


32.36 


13.2 


52.5 


98.0 


93±10 


13S 


1.824 


57.42 


23.9 


73.9 


155 


141±16 


14S 


2.278 


87.07 


37.4 


91.8 


216 


204±23 


15S 


2.795 


118.3 


52.6 


104 


275 


268±24 


16S 


3.376 


148.8 


68.6 


111 


329 


367±33 


17S 


4.024 


176.3 


85.0 


114 


375 


499±55 


IBS 


4.731 


200.9 


101 


113 


415 


614±59 


19S 


5.518 


221.5 


117 


109 


448 


741±66 


20S 


6.378 


238.1 


133 


104 


475 


1007±110 




8D 


0.474 


12.70 


4.19 


5.7 


22.6 


16± 7 


9D 


0.665 


52.76 


17.4 


15.1 


85.3 


73±14 


lOD 


0.909 


138.0 


45.6 


28.4 


212.0 


188±25 


IID 


1.189 


268.8 


89.4 


43.8 


401.9 


369±38 


12D 


1.524 


424.3 


145 


58.0 


627.2 


633±67 


13D 


1.905 


603.4 


207 


69.8 


879.9 


830±69 


14D 


2.341 


770.4 


270 


78.3 


1119 


1125±128 


15D 


2.835 


920.0 


331 


83.4 


1335 


1309±70 


16D 


3.387 


1046 


389 


85.4 


1521 


1420±82 


17D 


4.001 


1147 


443 


85.1 


1675 


1630±67 


18D 


4.680 


1225 


492 


82.9 


1800 


1804±88 


19D 


5.423 


1280 


538 


79.3 


1897 


1978±96 


20D 


6.235 


1317 


580 


74.9 


1972 


2036±60 



n'P states (figure 5(a)). These states have significantly larger Wbbr values than that 
of the 16S state, so that population transfer to the n'P states may affect the resultant 
BBR photoionization signal. As an example, the calculated rates of spontaneous and 
BBR-induced transitions from the 16S and the 16D states to the n'P states are shown 
in figures 5(b) and 5(c). 

Another feature is that BBR also induces transitions to high Rydberg states 
(figure 5(a)), which are field ionized by the second electric field pulse extracting the 
ions to channeltron for detection. This 100 V/cm pulse ionizes Rydberg states with n > 
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Figure 5. (a) Schematic illustration of BBR induced and field ionization processes 
upon excitation of the 16S state. They include redistribution of population over other 
n'L' Rydberg states due to spontaneous and BBR induced transitions from the initial 
state, (b) Calculated spontaneous and BBR induced transition rates from the initial 
16S state to other n'P states, (c) Calculated spontaneous and BBR induced transition 
rates from the initial 16D state to other n'P states. 



40, according to equation (jHl). Numerical calculations showed that high lying Rydberg 
states are populated at rates comparable with the direct photoionization rate Wbbr 
given in Table 1. Since all the states above Uc = 40 are field ionized, all BBR populated 
states with n' >40 will contribute to the detected Na"*" ion signal. The total rate Wr 
of BBR transitions to all Rydberg states with n' > 40 was calculated by summing the 
individual riL n'L' contributions given by equation (jH)): 



R 



J2 Yl WinL^n'L' 



(15) 



n'>40 L'=L±l 
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The calculated W^i values are given in Table 1. Note, that these ionization rates 
amount to between 30 to 50 % of the direct photoionization rate values, Wbbr, and 
they obviously cannot be ignored. 

Another kind of processes that may also lead to the redistribution of populations 
between Rydberg states is uL mixing in collisions with the abundant ground state atoms. 
However, the estimated rates of this process for a beam of density n^s ~ 2 x 10^'^ cm~^ 
produced by a source at 635 K were found to be much lower than the rates of spontaneous 
and BBR-induced transitions. Therefore the contribution of collisional uL mixing can 
be safely disregarded. 



4-3. Total BBR ionization rates 

When the rates of redistribution of population over n'V Rydberg states following laser 
excitation of a single nL Rydberg level are known, the total ionization rate leading to 
the production of Na"*" ions can be calculated. Assume that some nS state is initially 
laser excited. If the rate of Penning ionization is negligible, the production of Na"*" ions 
is described by the rate equation 

'^-^^ = Wf^nNns it) + E W-'i^N^^p it). (16) 

n' 

Here, IV^Bi? — ^bbr + Wr is the total BBR ionization rate of the laser excited nS state, 
and it consists of two contributions: the direct BBR photoionization rate Wbbr of the 
nS state given by equation (fTT|l . and the rate Wr of equation (|T5|l for BBR induced 
excitation of states with n > 40, which are subsequently field ionized by the second 
extraction field pulse. The main contributions to the sum in equation ()16p are those of 
n'P states with n' = n, {n — 1) , which are the nearest neighbors to the initial nS state 
(see figures 5(b) and 5(c)). The time evolution of the number of n'P atoms is described 
as follows: 

^^^^ = W{nS^ n'P) N^s (t) • (17) 
at 

Solving equations (fT^ and (fTTj) . we obtained the total effective ionization rate of the 
Na(nS) atoms: 

K%R = Wbbr + Wn + W^™., (18) 

where Wmix is a contribution from the spontaneous and BBR-induced mixing of 
neighboring states: 

't2-ti 



Wmix = Y,W{nS^ n'P) WiE^ 



1 ^e//. (19) 



Here te// is the effective interaction time given by equation (jSJ, and Te// is the effective 
lifetime of the laser excited nS state. Similar equations were derived for nD states, 
except that Wmix accounted for the transitions to both n'P and ra'F states. Results of 
numerical calculations using the semi-classical formulae for matrix elements from 
are summarized in Table 1. 
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Effective quantum number 

Figure 6. Total BBR induced ionization rates for (a) nS states and (b) nD states. 
Open circles and squares - experiment, solid lines - theory. 

4-4- Comparison of theory and experiment for BBR ionization rates 

Experimental and theoretical data on the total BBR photoionization rates are compared 
in figure 6. The solid curves are the theoretical values of W^^j^ from Table 1. The 
squares and circles are the corresponding experimental values M/^^rr obtained using 
equations (jH) and ((Zj) (in equation (jlj) the value of Wbbr was replaced with 
Experimental data were averaged over 5 measurements. The normalization coefficient 
C/) in equation ((Zj) was the only parameter whose absolute value was adjusted to fit the 
experiment to the theory. A good agreement between the experimental and theoretical 
data was found for nD states (figure 6(b)). However, the data for nS states, obtained 
with the ratio C £,/ C 5=1.6 measured earlier, exhibit considerable discrepancies for states 
with higher n (figure 6(a)), while a better agreement for states with lower n is seen. 
The values of W^^rr exceed the values of W^^br by 2.1 times for n=20, and the shape 
of the experimental n dependence is significantly different from the theoretical one. 

One possible explanation of such anomaly for nS states is related to their specific 
orbit that penetrates into the atomic core. The penetration causes a strong interaction 
between the Rydberg electron and the core, e.g., due to core polarization jHSj- This 
results in a large quantum defect and in a Cooper minimum in the photoionization cross- 
sections. Each nS state has its individual photoionization cross-section and position of 
the Cooper minimum. The accuracy of theoretical frequency-dependent cross-section 
values affects the total BBR ionization rate when integration over the Planck distribution 
in equation (fTTjl is performed. It is likely that the accuracy of theoretical models for nS 
states is affected by difficulties to account for the core penetration. This assumption is 
supported by good agreement of theory and experiment for the hydrogen-like nD states, 
which have a small quantum defect and an almost non-penetrating orbit. 

We conclude that our theoretical calculations of W^ni.ri are correct to better than 
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Figure 7. Experimental ratio R of Na^ and Na+ ion signals as a function of the 
effective quantum number n^ff for (a) nS states and (b) nD states. 



10% for nD states and lower nS states, and probably incorrect for higher nS states. 
New more accurate calculations of S state ionization rates are required if equation (jlOj) 
is to be used for the determination of absolute AI rate constants k^j for nS states. 
However, such calculations require development of new approaches that are not yet 
available for Rydberg states. Until new theoretical data become available, we rely on 
our experimental data for nS states instead of theoretical BBR ionization rates. The 
experimental values of W^^p^ are summarized in the last column of Table 1. These are 
our recommended values of the BBR ionization rates at 300 K, which should be used in 
equation (fTIH) . 

5. Experimental results for associative ionization rate constants 
5.1. Ratio of molecular and atomic ion signals 

Figure 7 shows the ratio R of equation ^ averaged over eight independent 
measurements at the density of atoms in the atomic beam 77,35=2 x 10^° cm~^. For 
the lowest states, the values of R reach 2.5 and 1.4 for nS and nD states, respectively. 
In both cases, R drops off rapidly and asymptotically approaches the value of 0.06 for 
highest n. 

Such nearly constant n dependence of R at higher n was unexpected. In order to 
verify whether the origin of the measured signal is collisional, the dependence of R on 
the atomic beam density was measured for low (9D and lOS) and high (19S and 18D) 
Rydberg states (figure 8). For the 9D and lOS states, R has a hnear dependence on 
nss, and hence the Na^ signal has a quadratic one. For 18D and 19S states, R is almost 
independent of n^Si hence the dependence of Na^ is linear. This means that for high 
n the apparent Na^ signal was not due to the associative ionization. Instead, it was 
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Atomic beam density (10 cm") 



Figure 8. Experimental dependences of the ratio R of Na^ and Na+ ion signals on 
the density n^s of ground state atoms in the beam for the (a) lOS state, (b) 9D state, 
(c) 19S state, and (d) 18D state. 



proportional to Na+ signal, and should be therefore considered as some kind of technical 
noise originating from Na^ ions. 

Analysis of the experimental arrangement used in the present study has shown that 
the noise was related to the Na"*" ions produced during the second extraction electric 
field pulse (figure 3(b)). The pulse acted on these ions for a shorter time than on 
Na^ ions produced before the electric pulse. Therefore, the former ions reached the 
channeltron at a later time than the latter ones. As a result, the detection gate, which 
was set for registration of the molecular ions, counted also the Na^ ions produced during 
the second extraction field pulse. This noise could not be eliminated by any technical 
means, and unavoidably had to be subtracted from the total ion signal counted during 
the molecular ion registration gate, although it was never reported in earlier experiments. 
The magnitude of this noise signal can be expressed by the formula 

Nnoise = N^L (0) P^AbR C-'^'^^ffAt, (20) 

where At is the time interval of the second electric extraction pulse, during which those 
Na^ ions are created, which can fall into the gate for molecular ions, and the time t2 
corresponds to the beginning of the extraction pulse. The noise is proportional to Wj^^j^, 
which grows rapidly with increasing n (figure 6). The time interval At was calculated 
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from the geometry of the detection system and the duration of the gate counting the 
molecular ions, and it was found to be about 150 ns. To find the contribution a of the 
noise to the value of R, equation pUjl must be divided by equation (jH) for Na+ ions: 

When this contribution is known, the absolute rate constant for AI in 
Na(nZ)+Na(3S) collisions can be extracted from the experimental ratios R of the 
molecular and atomic ion signals using equation (fTUI) modified by the correction a of 
equation (j2H): 

■u/tot 

kAi = {R-a)-^^. (22) 

The values of a for nS and nD states were calculated using equation (j2ip . 
For highest n of the experiment they practically coincide with the measured ratios 
(figure 8(c) and 8(d)) 



5.2. AI rate constants 

The experimental dependences of kAi on ne// were obtained from equation using the 
measured ratios of Na^ and Na"*" signals and the data of Table 1. These dependences are 
shown in figure 9. For the sake of completeness, we present two sets of results. Figure 
9(a) shows the rate constants determined with the experimental values W^^j^ of the 
BBR ionization rate, while figure 9(b) shows the same but using the theoretical W^^j^ 
values. The plots in both figures are nearly identical, except for some minor differences 
in the absolute values of kAi- Noteworthy, all dependences show the maximum in the 
AI rate near n^ff ~ (11 — 12), and fall off rapidly with increasing ne//. 

Comparing our results with those of earlier AI studies, shown in figure 1, we see 
a noticeable difference with the earlier data of ^2] on k^i for nS and nD states which 
were also obtained in a single beam experiment: (i) in our study the rate constants of AI 
for nD states are nearly twice larger than those for nS states, while opposite ratio was 
reported in ^2]; and (ii) our values of kAi fall off more rapidly with increasing n than 
those of ^21- We believe that our measurements are more reliable, since, in contrast to 
[12] , using equation we did not need to measure the number of Rydberg atoms and 
the effective lifetimes, and thus we could avoid the related experimental uncertainties. 

The only important parameter, which was not exactly known in the present 
experiment, was the density n^s of ground state atoms in the beam, which was calculated 
using the vapor pressure formula of Nesmeyanov |21]. At the same time, care was taken 
to ensure that the source temperature and the density of atoms in the beam were 
constant during the experiment. Therefore, the relative values of kAi for different nS 
and nD states are not affected by this uncertainty. 

Another interesting observation can be made by comparing our data of figure 9(b) 
with the n dependence for nP states of figure 1. Such comparison is given in figure 10. 
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Figure 9. Experimental data on the associative ionization rates for nS and nD states, 
obtained from equation (|22() using the measured ratio of Na^ and Na"'" ion signals, and 
(a) the experimental values of W^^j^ ; (b) the theoretical values of Wggj^ 



One can see that both data sets have similar absolute values and similar dependences on 
the effective quantum number ^e//. Note, however, that the data reported for nP states 
in ^21 consist actually of the AI rate constants for the states with 5<n<15, which were 
measured in jH] under quite different conditions of a crossed beam experiment, and the 
rate constants for n > 15, which were actually measured in a single beam experiment 
. Joining of the two data sets was used in ^2] for scaling their relative cross section 
values for n > 15. 



6. Discussion 



Our measurements of associative ionization rate constants were based on use of 
equation (j22I)- The accuracy of the measured kAi values depends on the accuracy 
with which the values of parameters R, a, W^^^, and n^s are known. The ratio R 
of Nag" and Na+ ions was measured directly and was found to be independent of laser 
intensities, detunings, and effective lifetimes. The measured n dependences of figure 7 
were reproduced in several experiments, which allow us to be confident that the values 
of R are determined correctly. However, with increasing n the fluctuations of R became 
important, because equation uses a product of (R-a), which decreases for high 
n to values close to 0, and W^^j^, which rapidly grows with increasing n. Therefore 
the reported values of kAi can be considered as reliable for lower n<16, whereas for 
higher n the error of determination of k^i is relatively large. More reliable data for 
higher n could be obtained by increasing the density of atoms in the beam by an order 
of magnitude or more. Unfortunately such increase inevitably leads to complications 
in the analysis due to additional phenomena, like radiation trapping, which was found 
to strongly affect the velocity distributions of excited atoms in a beam experiment at 
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Figure 10. Comparison of tlic experimental associative ionization rate constants 
for nS (open circles) and nD (open squares) states obtained in the present study 
(figure 9(b)) with the experimental rate constants for nP states (solid circles) obtained 
in and [El (figure 1). 



densities of 2x10^^ cm"^ ^7] . 

Another possible source of errors in the determination of kAi is the choice of values 
of the total BBR induced ionization rate W^^j^. Figure 6 shows a good agreement of 
experimental and theoretical values for all nD states, while for higher nS states the 
agreement is rather poor. Therefore, one must choose whether the theoretical or the 
experimental values of BBR ionization rates, given by the last two columns of Table 1, 
will be used in equation (j22I)- Since the experimental results of figure 6 were reproduced 
in many experiments, we believe that the experimental values W^^^ are more reliable. 
Therefore, figure 9(a) should be considered as the final result of the k^j measurements. 

The absolute values of k^i depend also on the absolute value of n^s calculated 
using the Nesmeyanov's formula [25 ■ Numerous studies involving alkali atoms have 
confirmed its accuracy in predicting the vapor pressures of alkali gases (see, e.g., the 
argumentation in [SH]) within a few tens of per cent. Note, that the relative values of 
kAi were not affected by the absolute value of n^s, since it was kept constant during 
the experiment. We conclude that on average our absolute values of kAi are accurate to 
within better than 50%, while the relative values of k^i are accurate to within (10-20)%. 
However, due to the specific n dependence of the term {R — a) in equation (j221), the 
reliability of k^i values is better for lower n. 

In conclusion, we have measured, for the first time, AI rate constants for sodium 
nS and nD states in the range of 8<n<20. In this study, we have implemented a 
novel method, in which the absolute values of rate constants were determined from the 
measurements of relative ratios of the molecular and atomic ion signals, with atomic 
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ions being from the BBR ionization of Rydberg atoms. The main advantage of this 
approach is that the measurements of absolute populations and effective lifetimes, which 
are associated with large experimental uncertainties, can be avoided. 

To implement this method, we have measured BBR ionization rates for nS and nD 
states of sodium with n = 8-20. Comparison of these data with the theory showed a good 
agreement for nD states, and disagreement for nS states with n > 15. This discrepancy 
cannot be attributed to any of the experimental uncertainties, and the theory must 
substantially revise the calculations for nS states. 

The comparison of our experimental data on AI rates with other experiments for 
nS and nD states is hampered by lacking overlap of the studied ranges of n. However, 
for the high n range, where some overlap with the data of [12] exists, notable differences 
were observed. This may be explained by the way we accounted for the background 
Na^ ions. These ions were created due to ionization during the extraction electric field 
pulse, and they arrived to the detector simultaneously with Naj . Such possibility was 
not considered in earlier studies [01 CHI HH IE], and it is not possible to reconstruct 
whether such background was present in these experiments. 

Curiously, our results for nS and nD states are in a good agreement with the data 
of [H] for nP states, though that experiment was performed with two crossed atomic 
beams. In the present single-beam study, ionization is due to head-tail overtaking 
collisions between ground-state and Rydberg atoms. Accordingly, the collision energies 
are about three times lower than those in the crossed-beam experiment, which used 
beam sources of the same temperature. In the second part of our article we provide 
a detailed comparison of the existing experimental data on associative ionization rate 
constants with the theoretical calculations for various relative velocity distributions. The 
calculations are based on the stochastic ionization model [H El El IHl EI , and they take 
into account all possible collisional, spontaneous, and BBR-induced mixing processes. 
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